Few emitters in a cavity: From cooperative emission to individualization 
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We study the temporal correlations of the field emitted by an electromagnetic resonator coupled 
to a mesoscopic number of two-level emitters that are incoherently pumped by a weak external drive. 
We solve the master equation of the system for increasing number of emitters and as a function of 
the cavity quality factor, and we identify three main regimes characterized by well distinguished 
statistical properties of the emitted radiation. For small cavity decay rate, the emission events are 
uncorrelated and the number of photons in the emitted field becomes larger than one, resembling 
the build-up of a laser field inside the cavity. At intermediate decay rates (as compared to the 
emitter-cavity coupling) and for few emitters, the statistics of the emitted radiation is bunched and 
strikingly dependent on the parity of the number of emitters. The latter property is related to the 
cooperativity of the emitters mediated by their coupling to the cavity mode, and its connection 
with steady state subradiance is discussed. Finally, in the bad cavity regime the typical situation of 
emission from a collection of individual emitters is recovered. We also analyze how the cooperative 
behavior evolves as a function of pure dephasing, which allows to recover the case of a classical 
source made of an ensemble of independent emitters, similar to what is obtained for a very leaky 
cavity. State-of-art techniques of Q-switch of resonant cavities, allied with the recent capability to 
tune single emitters in and out of resonance, suggest this system as a versatile source of different 
quantum states of light. 

PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Gy, 42.65.Hw 



I. INTRODUCTION 

The interaction of an ensemble of identical two-level 
emitters with the electromagnetic field has been the sub- 
ject of intense theoretical and experimental studies since 
the pioneering work of Dicke Among the others, su- 
perradiance (or superfluorescence) is the most celebrated 
effect; it is manifested as a transient effect by the delayed 
emission of an intense burst of light, after the atomic 
population has been totally inverted [5]. In essence, su- 
perradiance is characterized by the emergence of collec- 
tive coherence, which leads to a strongly enhanced emis- 
sion rate as compared to the case of uncorrelated emit- 
ters. Such striking consequences of the quantum interfer- 
ence effect in radiation-matter coupled systems has been 
the subject of an intense experimental interest both for 
atomic [21 E] and solid state [4] ensembles. Early stud- 
ies on the statistical properties of the emitted light have 
evidenced the role of quantum fluctuations in the initial 
stages of the superradiant emission [MZI • Whereas super- 
radiance usually appears in the transient regime, recent 
theoretical studies have evidenced other signatures of col- 
lective behavior in the steady state regime, e.g. when the 
two-level emitters are continuously excited by incoherent 
light [8HTT]. In these references, cooperative emission is 
obtained by coupling the emitters to a low quality (Q) 
cavity mode, which behaves as a collective channel of 
losses. As a matter of fact, strong confinement enhances 
the interaction of each emitter with the cavity field, even- 
tually making it predominant with respect to individual 



relaxation channels. Collective dynamics can be induced 
this way, provided the atoms identically interact with 
the cavity mode, a condition much less stringent than 
in free space. Tuning the pumping power allows one to 
explore different regimes, leading to different intensities 
and statistics for the emitted radiation [101 [H]. In partic- 
ular, it has been shown that when the pumping strength 
is too low to saturate a single emitter of the ensemble, 
the collection of two-level systems emits a highly bunched 
field f^llll]. The strong temporal correlation has been at- 
tributed to the relaxation of the system in the subradiant 
states of the ensemble of emitters. These states behave 
as sources of photon pairs triggered by the driving pulses, 
hence the bunching effect in the correlation properties. 

In this paper, we analyze the subradiant regime and 
go beyond previous works by keeping the cavity mode as 
a relevant degree of freedom. This situation is new with 
respect to previous theoretical studies, where the radi- 
ation mode was commonly eliminated by an adiabatic 
approximation. We show that the bunching degree de- 
cays with the number of emitters and oscillates with its 
parity, being larger when the number of emitters is even, 
upon converging to a single value as the number of emit- 
ters grows. We give an interpretation of this effect, based 
on the structure of Hilbert space of the system. Further- 
more, keeping the cavity mode in the description allows 
to explore two other limits, depending on the cavity fi- 
nesse. If this parameter is significantly increased, the 
strong coupling conditions lead to the emission of a field 
showing no temporal correlations, an effect that we iden- 
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FIG. 1: Schematic of the system under consideration in this 
work: TV two-level emitters incoherently pumped by an exter- 
nal drive, coupled to a single resonator mode. Correlations 
in the emitted field can be detected through a HBT set-up 
either from the cavity or from the emitters' loss channel, re- 
spectively. A typical incoherent pumping scheme is shown, 
in which high-energy excitations relax down to a well defined 
ground state transition (e.g. through phonon scattering) via 
photon emission. 



tify with the emission of Poissonian radiation similar to a 
lasing transition. Recently, such regime has been investi- 
gated from the spectral point of view [12 [T3] , while here 
we address its statistical properties. On the contrary, for 
decreasing cavity finesse one progressively decreases the 
coupling with the mode, allowing to continuously explore 
the boundary between collective and individual behavior 
of such radiation sources. When all the emitters become 
uncorrelated, increasing their number one by one up to 
large values coincides essentially with building a classical 
source with chaotic statistical properties. Dipole corre- 
lation can also be destroyed with pure dephasing, a situ- 
ation that we eventually explore as well. 

The present work becomes particularly relevant when 
considering that the subradiant regime is reachable if 
and only if each single emitter of the ensemble is pre- 
dominantly coupled to the cavity mode. In particular, 
this allows to pump each emitter very slowly, but still 
faster than any individual relaxation mechanism. This 
is nowadays possible thanks to impressive progress in 
nano-fabrication technology on single emitter-cavity cou- 
pling [Tm27j . Small assemblies of artificial quantum 
emitters, like quantum dots |14j and superconducting 
qubits |15j coupled to a high-Q resonator mode have re- 
cently become available, strongly motivating to investi- 
gate the new physics provided by these systems. More- 
over, beyond introducing an external parameter to con- 
trol the steady state collective regime, this model is of 
great interest also for the study of lasing transition, as 
well as to investigate the connections between superradi- 



ance and lasing [25]. 

The paper is organized as follows. We first introduce 
the model and the methods employed to solve it, both nu- 
merically and analytically (Sec. II). Then, we present re- 
sults as a function of cavity finesse, and we give a physical 
interpretation for the different regimes found (Sees. Ill, 
IV, and V). Before summarizing the main conclusions, 
we give a brief description of possible realistic implemen- 
tations to explore the physics investigated in this work 
(Sec. VI). 



II. MODEL AND METHODS 

The system under study consists of iV two-level emit- 
ters coupled to the mode of an electromagnetic cavity of 
annihilation (creation) operator a (a^) with equal cou- 
pling strength, g, and it is schematically pictured in 
Fig. [l] To focus on collective effects, we will assume 
that the resonator and all the emitters have the same 
transition frequency, w. The range of validity of this 
ideal model in practical implementations is discussed in 
Sec. VI. The model can be rigorously described by the 
Tavis-Cummings Hamiltonian (in the rotating wave ap- 
proximation) 

N 

H = '^huj{al) + hoja^a + ig{aJ+ - J^a^) , (1) 

i=l 

where J+ — cr!|_ is the collective dipole operator, and 
a'_ = |gj)(ej| is the lowering operator of the single emit- 
ter i, whose excited (ground) state is denoted by je^) 
{\gi))- The Hilbert space corresponding to the assem- 
bly of resonant atoms is equivalent to the tensor prod- 
uct of the 2-dimensional space of each single emitter, i.e. 
v. ~ (C^)®^, spanned by the basis of factorized states in 
the form |ii) |«2)-"|*n)) where = {5n,e„} (in the follow- 
ing we shall refer to this basis as the uncoupled basis). 
When including possible dissipation channels and exter- 
nal driving, the Markovian dynamics of the full system 
is given by 

N 

p = ^^[H,p]+Y,^^ (P) + (p) + Cy- ip) ■ (2) 
1=1 

The Lindblad dissipative terms, modeling incoherent loss 
and pumping of each emitter, explicitly read 

A(p) = C^{p) + Cp^ip) 

- ^(2fTLp<-<aLp-pa>L)-f 

^(2a>aL-aL<p-paK), (3) 

where the quantity 7 is the individual emission rate in 
loss channels distinct from the cavity mode for each emit- 
ter, and Px quantifies the excitation rate. As pictured 
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in Fig. [Tj incoherent pumping can be obtained by reso- 
nantly exciting higher energy states, that further relax 
down to the levels of interest. The decay of the cavity 
mode with rate k through imperfect mirrors is described 
by the following term 

^a{p) = ^(2apa"f - a^ap- pa' a) (4) 

whereas the pure dephasing of the emitters, which is un- 
avoidable in solid-state artificial atoms, can be taken into 
account by a Liouvillian term of the form '29' 

i=l 

Following standard input-output theory [30], the out- 
put field at the photodetectors is proportional to the a 
operator (and its conjugate). We are interested in the 
temporal correlations of the field radiated by the cav- 
ity mode. Experimentally, an ordinary Hanbury-Brown 
Twiss (HBT) setup can be used, giving access to the 
steady state second-order correlation at zero-time delay 

m 

5i'^(0)-.gi?^-(aVaa)/(ata)2. (6) 

In order to calculate the latter quantity, we numerically 
solve Eq. ^ either by direct integration or by Quantum 
Monte Carlo (QMC) simulations Direct integra- 

tion allows one to obtain accurate results at lower com- 
putational time costs, for which we use an extension of 
the code already employed in a previous work [34' . How- 
ever, the tensor structure of the phase space makes the 
computational resources needed for a direct integration 
rapidly diverge. On the contrary (wave function) QMC 
gives the solution as a stochastic average over quantum 
trajectories of wavefunctions unravelings, drastically re- 
ducing the amount of space needed and allowing us to 
push forward the previous computational limits |33j . 

If K is large enough, the cavity mode can be adiabat- 
ically eliminated from the equations, giving rise to an 
effective dynamics described by the reduced master equa- 
tion 

N 

p = £r(p) + EA(p), (7) 

i=l 

where 

Cr{p) = ^{2J-pJ+-J+J-p-pJ+J-). (8) 

The physical meaning of this equation is the following: 
the cavity acts as a common relaxation channel for the 
atoms, whose dissipation rate is quantified by the cou- 
pling strength with the collective atomic dipole J_, i.e. 
r = Ag^ / K. Namely, adiabatic elimination is valid when 
K ^ NT [10]. In this regime, the cavity field, and thus 
the output field is proportional to the collective dipole 



degree of freedom, J„. Consequently, the temporal au- 
tocorrelation function at zero time delay can be expressed 
as 

fff (0)^(J+J+J_J_)/(J+J_)2. (9) 

Note that the same detection properties could be ob- 
tained by registering the far field where the atomic 
dipoles interfere, as it is pictured in Fig. [l] Also for 
Eq. Q we can employ numerical approaches either 
based on direct integration or on QMC, respectively. As 
usual, once the master equation is solved for p the ex- 
pectation values can be calculated with the general rule 
(O) — Tr{Op}, where O identifies the operator corre- 
sponding to a given observable. 

Sweeping the cavity finesse allows to continuously ex- 
plore the transition from a collective behavior driven 
by Ct{p) to an individual one, corresponding to Ci{p). 
In this work we analytically study some limits of these 
regimes. If the dynamics is individual, it is naturally 
described in the uncoupled basis mentioned above. If 
the dynamics is collective, an order appears in the sys- 
tem due to the build up of correlations between dipoles 
belonging to different atoms [5], manifested by the oc- 
currence of coherences in the uncoupled basis. In the 
latter natural basis is provided by the eigenstates 

|J, M) of the total angular momentum {J^, J^} (in the 
following, we shall talk about the coupled basis of the 
atomic states). Note that among them, the states satisfy- 
ing J_ I J, M) = are not coupled to the electromagnetic 
field and are thus defined as "subradiant" . 



III. SECOND-ORDER CORRELATIONS OF 
THE EMITTED FIELD 

As mentioned above, the regime studied in the liter- 
ature so far is the adiabatic elimination regime, where 
the cavity mode operator can be assumed as a oc J_. 
We assume a fixed and weak incoherent pumping rate, 
Pxl 9- In order to understand the collective behavior and 
to link it to the previous literature, we plot in Fig. [2^ the 
average cavity mode occupation, Ua = {a^a), the average 
collective occupation of the emitters, nj = (J^J^), and 
the total emitters population, Nux (where = {a\_a^_) 
in the steady state. The second-order correlations are 
plotted in Fig. [2]d, and they have been calculated both 
from the cavity mode operator, Eq. ([6|, and the collec- 
tive dipole operator, Eq. respectively. The results 
are shown for N = 5 emitters in the cavity, which repre- 
sents a good compromise to simultaneously have a well 
converged numerical solution for Eq. Q and to evidence 
the effects of collective behavior. As it can be seen in Fig. 
[2]d, the two calculated correlations don't match for most 
of the parameter space, showing that the cavity mode 
has to be taken into account to obtain the correct result 
when the cavity finesse is large. The two curves converge 
to the same value when K/g ~ 50, where the adiabatic 
approximation starts to be satisfied. 
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FIG. 2: Full numerical results for average populations and 
second-order correlation in steady state as a function of k 
(normalized to g) for the case of A'^ = 5 atoms in the cavity, 
calculated from atomic and cavity fields, respectively. Param- 
eters are: = O.OOI5, 7 = O.OOOlp. 



We first focus on the case where adiabatic elimination 
is valid: two different situations emerge. The case where 
r = 4(7^/«; > 7 corresponds to a collective dynamical 
behavior for the atomic ensemble. In the low pumping 
regime under study, the emitted field is strongly bunched, 
(7^^-'(0) > 1, confirming previous theoretical works [SHU] 
on a similar model. In the two-atom case, this behavior 
has been attributed to the emission of so called "super- 
radiant" photon pairs [9], typically separated by a delay 
scaling like the inverse of the pumping rate. This regime 
has been further explored for a larger number of emitters 
in [ini E] , where the authors show that the atomic en- 
semble has mostly relaxed in the subradiant states. This 
is in agreement with Fig. [2^ where the average inten- 
sity emitted by uncorrelated atoms, Nux is compared 
to the average fluorescence field emitted by the collec- 
tive atomic dipole, {J+J-). It can be clearly observed 
that in a n/g range spanning several decades, the con- 
dition (J+J^) < Nux is usually satisfied, a signature of 
the anticorrelation of the atomic dipoles and of the ef- 
ficient relaxation in the sub-radiant states. Note that 
the subradiant regime is reached when < T, a condi- 
tion that is not fulfilled anymore for n/g > 10"^. In that 




10""* 10"^ 10"^ 10"^ 10° 10^ 10^ 10^ lO"* 10^ 10*^ 



K 

FIG. 3: Full numerical solution for zero time delay second- 
order correlation in steady state calculated from the cavity 
degree of freedom, as a function of k (normalized to g), for low 
incoherent pump rate Px = O.OOlgr and emitters' dissipation 
7 = O.OOOlf?, and varying number of emitters, A'^. Left inset: 
oscillations as a function of A'^ in the good cavity, subradiant 
regime. Right inset: bad cavity limit of correlations as a 
function of A^. 



case the steady state emission becomes super-radiant, i.e. 
(J-(-J_) > Nux, and correspondingly g'^''(0) approaches 
1, a limit that is exactly reached only for ^ 1 as 
evidenced in previous work Finally, when n/g is 

increased up to a point where 7 > F (very bad cavity 
regime) , emission appears as if it was coming from an in- 
dependent collection of emitters. Under such conditions, 
the limiting values are such that (J+J_) — > Nrix and 
5(^^(0) 2(1 — 1/iV), a limit that we analyze in detail 
below. 

We focus now on the statistical properties of the system 
by looking at the dependence on the number of emitters, 
the results are shown in Fig. [3] with the same parameters 
used in Fig. [2] Keeping the cavity mode as a relevant 
degree of freedom allows one to explore the temporal cor- 
relations on the entire n/ g range, which can only be done 
by the full numerical solution of the master equation. At 
this point, we go beyond previous works. First, we no- 
tice that a transition to amplified spontaneous emission 
clearly occurs for very large cavity finesse, i.e. as soon as 
n/g < 0.001 {Px > k), with the temporal fluctuations be- 
coming Poissonian independently of the number of emit- 
ters. A temporal picture for this effect is that as soon 
as Px > K, the cavity stores the photons on a timescale 
longer that the typical delay between two emission events 
by the atomic assembly, blurring out any temporal corre- 
lation feature. This situation has been already evidenced 
for the single emitter case in [5^36j . and it was identi- 
fied as a lasing threshold signature. Here, the system 
becomes a source of radiation with Poissonian statisti- 
cal noise, reminiscent of a few emitters laser. The case 
for N = 2 in the strong emitter-cavity coupling regime 
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was also studied for a similar model in Ref. [371 131] as a 
function of pump power. As a matter of fact, three main 
effects on measurable quantities can be identified: (i) the 
average cavity population becomes larger than one, set- 
ting the onset of stimulated emission (see also Fig. [2^) , 
(ii) the population of each emitter is clamped to a def- 
inite value, here, ~ 0.5 (i.e. Sn^, — >■ 2.5 in Fig.^k), 
and (iii) g'^HO) ~^ 1 (more clearly evidenced in Fig. ^ffor 
any N = 2,3,4,5). 

On increasing n/g, it is evident from Fig. [s] that in 
the subradiant range there is an unexpected and quite 
pronounced parity-dependent oscillation of the photon 
correlation value at zero-time delay, a behavior that has 
not been discussed in the literature, to our knowledge. 
Such parity oscillations are not limited to the adiabatic 
elimination region (see the inset on the left panel for 
K/g — 0.45 and up to a number of emitters = 8). 
The latter calculation required consistent numerical op- 
timization and was performed on a cluster, where the full 
master Eq. ^ could be solved without truncation effects 
on the convergence. Although quite difficult to be ana- 
lyzed in the dressed state picture, the oscillatory behavior 
of correlations as a function of the number of emitters can 
be roughly interpreted as a consequence of the collective 
coupling in Eq. ([s]), which dominates the dynamics over 
the individual relaxation channels in Eqs. ([s]) and Q. 
The former couples the states with the same quantum 
number J and is diagonal in the coupled basis. On the 
contrary, the other channels act on each emitter and are 
responsible for coherences between states with different 
values of J. We will investigate this regime in detail in 
Sec. IV. The bunching behavior reaches a maximum for 
K ~ slightly depending on the number of emitters. An 
interesting even/odd behavior has also been found in the 
spectral response [12] and it would be interesting to re- 
late the two phenomena in future developments of the 
present work. 

Finally, in the opposite limit of n/g becoming large 
enough that 7 > F the radiation is characterized by a 
field emitted from independent atoms. Oscillations tend 
to disappear, as it can be seen in Fig. [s] (see also inset 
on the right hand side). In this regime, the single emit- 
ters undergo individual dynamics, and the dipoles be- 
come uncorrelated (see also the large k limit in Fig. [2^, 
where {J^J^) — > Nrij.. Restricting our considerations 
to the correlation function, an analytical expression of 
5^(0) is derived in Sec. V, where we show that ^^^•'(O) = 
2(1 — 1/N). Note that in the limit of a large number of 
atoms, 5^^'(0) — 2, which corresponds to the emission of 
an incoherent chaotic field. 



IV. PARITY-RELATED OSCILLATIONS 

In this Section we calculate analytical expressions for 
the bunching degree of the emitted field, in the two and 
three emitters cases, which provide a physical explana- 
tion for the observed parity dependence. The analytic 



treatment is valid in the regime of adiabatic elimination 
of the cavity mode dynamics. To sustain our analysis, 
we have also used a Quantum Monte Carlo (QMC) sim- 
ulation [531. We have calculated 5*-^^(0) as a function of 
the number of emitters N, with N being as large as 10. 
The results are plotted in Fig. 4^ by simulating the re- 
duced master equation, Eq. (17). As it can be seen, we 
alternatively recovered the result shown in Fig. [3] at a 
specific value of n/g — 100, by generalizing to a larger 
number of emitters that would be unaccessible through 
direct numerical integration. Again, bunching and os- 
cillatory behaviors are clear signatures of this collective 
low pumping regime, which we are going to analytically 
describe below. 



A. Two emitters case 

In the following of this Section we use the coupled ba- 
sis for the emitters. Such basis of states is schematically 
indicated in Fig. |4]d. In particular, the two-emitter case 
is pictured in the first panel of Fig. [4]d, and we will fo- 
cus on this diagram here. In this specific case, the basis 
I J, M) is non-degenerate and the Hilbert space consists 
of the three symmetric (or Dicke) states |1, 1) = |ei,e2), 
|1,0) = (|ei,52) + l5i,e2))/y2, |1,-1) - 151,32), and the 
antisymmetric state |0,0) — (|ei,52) — \gi,S2))/V^ (i.e., 
the usual triplet and singlet states of a spin 1 system). 
The steady state (and g'^^^O) statistics) is the result of a 
competition between collective (specifically Cr) and in- 
dividual decoherence channels (i.e., Lp_^,L^, C.^*). The 
results shown in Fig. |4^ are obtained without pure de- 
phasing and with = IO7. In the steady state, low 
pumping equally populates the subradiant states |1,— 1) 
and |0, 0). As it was first discussed by Temnov and Wog- 
gon |9j , one can decompose the problem in two mutually 
exclusive cycles. In the regime under study, the F-related 
de-excitation events are much more probable than the ex- 
citation events driven by the pumping rate P^. Hence, 
as soon as an excitation is stored in a state where the 
two mechanisms are possible, it is more likely that the 
excitation will undergo a collective emission than an up- 
ward pump jump. The microscopic scenario will then 
involve one cycle (pictured in green) consisting of the 
excitation and relaxation of the symmetric state |1,0), 
leading to the emission of single photons. In the other 
cycle (pictured in red), the excitation is first stored in 
|0,0). Then, during the time evolution, the only pos- 
sibility is to excite the state |1,1), leading to the fast 
emission of two photons. The system prepared in the 
subradiant state |0,0) behaves thus like a determinis- 
tic source of photon pairs. Each pair is emitted on a 
timescale l/F, the delay between two pairs being typi- 
cally 1/Px- Consequently, the emitted field is expected 
to show bunching, ^^^•'(O) '-^ F/Pj,, an order of magnitude 
that we exactly recover below by analytically solving the 
approximate master equation. 

In fact, to go beyond this qualitative interpretation we 
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have computed the expression for 5^^) (0) in terms of the 
relevant elements in the steady state density matrix 



(2)rn^ - 



Pl,l 

[pi,i + Pi,o] 



(10) 



at first order in P^j/r. The diagonal matrix elements of 
the density matrix are written with the row (column) in- 
dex only, e.g. (J, M|p|J, M) = pj^m- In this two-emitter 
case the evolution equations for the diagonal elements 
do not involve off-diagonal terms and the solution is ob- 
tained by solving a system of rate equations governing 
the evolution of the atomic populations. With lengthy 
but straightforward calculations, we express Eq. ([7| on 
the basis above, and solve it for the steady state p — Q. 
We finally obtain an analytical estimation 



16^ 

9 kP^. 



4 r 

9^ ' 



(11) 



which for Px/g = lO^'^ and n/g = 100 gives g^^^(O) ~ 
17.8, a value that is very close to the numerical results 
(obtained from both direct numerical integration and 
QMC). 
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B. Three emitters case 

In the iV = 3 case, the simple picture above can- 
not be applied directly. The most important difference 
emerges immediately from the very nature of the Hilbert 
space. As pictured in Fig. |4]d, right panel above, the 
J = Jmin = 1/2 subspaces have a 2-level system struc- 
ture. As a consequence, the subradiant states corre- 
sponding to \Jmim —Mj^.^^) do not bchavc as determin- 
istic sources of photon pairs as in the N = 2 case, but 
eventually give rise to single-photon cycles when they are 
pumped (light blue arrows in Fig. jija). The presence of 
these extra loops spoils the bunching of the emitted field, 
and this gives a simple explanation for the abrupt loss of 
g'^-* (0) bunching when going from two to three emitters. 

To sustain this argument we have computed the 
second-order correlation 5*^^^ (0) in terms of the density 
matrix elements, which for the = 3 case reads 



12(^2.2 +P2 i) 



[3p3 3 -f- 4p3 1 + 3p3 _ 1 + pi l.-i + pi 1 

'- 2'2 2'2 2' 2 2'2' 2'2 ' 

(12) 

We followed the same notation as in Eq. ( 10 ) for the 
matrix elements. The subradiant states are written 
|l/2, -1/2;1> and |l/2,-l/2;2) where the second 
label is a shorthand notation related to the g = 2 ex- 
tra degeneracy (see App. [A|and the related figure). Af- 
ter writing the master equation in this basis, we then 
proceeded with a perturbative solution of the linear sys- 
tem for the steady-state matrix elements in powers of 
(Px/T). Differently from the N = 2 case, the equations 
of motion for the diagonal terms now involve off-diagonal 



FIG. 4: (a) Oscillations in the second-order correlation at zero 
time delay in the bad cavity, sub-radiant regime {k/ g = 100). 
(b) Schematic energy level diagram for N = 2 and A'^ = 3 in 
the total angular momentum basis, and their generalizations 
for even and odd number of emitters, respectively. For N — 3, 
g = 2 reminds the two-fold permutation degeneracy 



matrix elements, preventing one from describing the dy- 
namics with a simple set of rate equations. Nevertheless, 
within the same perturbative picture, we can consider 
these terms to be much smaller than the diagonal ones. 
Indeed, the coherences are brought by some residual in- 
dividualization of the Pj.-driven channel that, for the pa- 
rameters chosen, acts as a perturbation over the collective 
dynamics induced by the F events. The latter argument 
has been checked with the steady-state density matrix 
elements obtained numerically. As detailed in App. [B] 
the final result reads 



(0) 



14 F P ( P- 

+ 2.067 - 3.501— + 9.18 — 

75 F If 



-23.4 



Px 



O 



Px 



(13) 



We obtain the analytical value ^'^-'(O) ~ 9.46, which 
compares very well with the numerically calculated value 
~ 9.55 reported in Fig. |4ji from QMC simulations. No- 
tice that the population of the state ||, |) cannot be ne- 
glected and was also taken into account in the analytical 
calculation. 
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C. Generalization to A'^ > 3 emitters 

Considering the energy level diagrams in the generic 
even and odd cases, as depicted in Fig. [4]d, lower panels, 
one can always find the two patterns above (i.e., N — 2 
and iV = 3) as distinctive features of each scenario. For 
increasing number of emitters more subradiant states be- 
come available, among which the population arising from 
the Px reservoir will be shared (bold levels in the figure) . 
This spreading of population unavoidably leads to the 
decay of the bunching degree as N is increased, as ob- 
served in Figs. |3] and |4^, and confirmed by the general 
trend anticipated in 9 where a semi-classical approach 
was used. In the limit of large large numbers, such as 
N = 10000, the bunching degree of g^^^O) converged to 
2, as for classical chaotic sources. Note that in our case 
taking the limit of large numbers requires some caution, 
as one could exit the validity of the adiabatic elimination 
approximation, k ^ NT, and eventually enter strong 
collective light-matter coupling regime (see, e.g., recent 
works [H1I12]). For the parameters chosen in the simu- 
lation in Fig. [4ji, the condition for adiabatic elimination 
limits the maximum number of emitters to about 2000. 



V. INDIVIDUALIZATION 

We consider now the case where the cavity finesse is 
decreased down to a point where the coupling of each 
emitter to the mode, F, becomes of the same order of 
magnitude of the coupling to its individual decay channel, 
7. As we have seen in Figs, and [3j this leads to totally 
uncorrelated emission events from the atomic dipoles. In 
this regime, which we define as the "individualization," 
the parity-induced oscillations vanish and the bunching 
degree converges to g^^\0) = 2(1 — 1/iV) (see Fig. [s] in- 
set on the right), a behavior that we demonstrate below. 
Note that the same value for the parameter 5*-^^(0) has 
been computed in the strong pumping case [TT1 . As a 
matter of fact, a large pump power leads to total pop- 
ulation inversion, a state where the dipoles necessarily 
become uncorrelated. Moreover, it is interesting to show 
in the context of the present work that an alternative 
way of obtaining a source of uncorrelated dipoles is to 
increase the pure dephasing rate, 7*. The latter effect 
is shown in Fig. [5] where a full numerical solution based 
on direct integration has been employed for n/g = 100. 
It is evident that individualization is reached as soon as 
7* > F, washing out the oscillations and leading to the 
same limit for the second order correlation parameter, as 
clearly displayed in the inset. 

For the sake of completeness, we have com- 
puted the steady state second-order correlation 5^^^(r), 
for any delay r. Adiabatic elimination is valid 
and this quantity can be calculated as a func- 
tion of the atomic operator, namely g^^' (r) = 
(J+(0)J+(r)J_(r)J_(0))/(J+(0)J_(0))2. We make the 
assumption that the steady state of the ensemble can be 



factorized in the uncoupled basis, which is reasonable as 
each emitter behaves individually. Thus, developing (see 
App. [C]) the expression of J± with respect to the indi- 
vidual atomic operator (Tj_ leads to 



E^9i'\r)lf + j:^^,lMl+griTrgrir)} 



(14) 

where we have introduced the intensity per emit- 
ter, the first and the second order correlation 
functions for the quantum emitter i, as li = 
(c7^(0)aL(0)), g^'\r) = (a;(r)aL(0))//„ ^^(r) = 
(cr^(0)cr!^(T)a!_(r)crL(0))//^, respectively. As initially 
mentioned in seminal lectures 1431 and papers [44i . and 



as it clearly appears in Eq. (14 1, the second order cor- 



relation function involves the emission of pairs of atomic 
excitations. When t — 0, this expression simplifies to 



g(2)(0) = 2(l-l/7V) 



(15) 



confirming our numerical simulations (see right inset of 
Fig. [2]). The bunching at small delays is due to an in- 
terference resulting from the indistinguishability of the 
atoms in the emitting pairs [43]. Note that surprisingly, 
this behavior was never observed experimentally. As a 
matter of fact, the interference is washed out as soon as 
individual first order coherence, (;'^^(t), vanishes [3Hl44l- 
W7\ . This phenomenon takes place on the timescale of 
dephasing processes, leading to the following expression: 



g(2)(0)^l-l/iV 



(16) 



eventually corresponding to a particle-like behavior. In 
fact, let's consider N independent sources each emitting 
classical identical particles one-by-one. Such a model- 
source is characterized by a total emission intensity / = 
J2i h with {Ii{t)Ij{t + r)) = a i ^ j. In the stationary 
regime we have {Ii{t)) = {Ii{t + t)) = /q, since we sup- 
pose that the N sources are characterized by the same 
main current. The total second order correlation signal 
is therefore 



{Iit)I{t + r)) 
{I{t)){Iit + r)) 



which becomes 



iV2 



(17) 



(18) 



However, our N idealized sources are single particle emit- 
ters i.e. gl (0) = 0. This clearly mimics the quantum- 
like behavior of two-level systems. Therefore, we get 
Eq. (16 1. Interestingly, from basic textbook quantum 
optics we obtain for a pure \n) photon Fock state the ex- 
pression (?'-^-'(0) = 1 — 1/n. The result is therefore also 
intuitive from this point of view. 

The limit given in Eq. ( 16 ) was already observed in 



the pioneering antibunching experiment by Kimble, Da- 
genais and Mandel with fiying sodium atoms coherently 
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FIG. 5: Zero time delay second-order correlation as a function 
of pure dephasing rate, 7* (normalized to g), at fixed Px = 
O.OOlg and k = lOOg. The inset shows the behavior as a 
function of at large 7*. 



driven |31| , and it was theoretically interpreted in the fol- 
lowing years [511 IHl - HT] . In these references, the source 
of dephasing was mainly attributed to the finite size of 
the atomic source with a random number (i.e. a Poisson 
distribution) of emitters inducing a spatial averaging out 
beyond the typical coherence volume of the fluorescence 
light [311 [46] . In nowadays experiments with solid state 
emitters confined in nanometre-size systems, this prob- 
lem could be ultimately overcome by considering a good 
control of the optical detection. However, notice that de- 
phasing processes are always present in solid state quan- 
tum emitters (due for example to coupling with phonons 
in solid matrices), and that the measurement of g^'^^O) 
involves actually an average over finite values of r near 
T ~ 0. If the typical dephasing time r* = I/7* is larger 
than the time resolution of the HBT setup (e.g. the jit- 
ter time of single photon detectors) this may explain why 
experiments reported so far with different quantum emit- 
ters H5H5^ tend to evidence a behavior that is closer 



to Eq. ( 16 ) rather than the present prediction described 
by Eq. (15 1. Observing the 2(1 — 1/A^) behavior, and 
its transition to 1 — 1/A^, could be a challenging issue, 
e.g. requiring a control over the local temperature of the 
emitters. 



VI. PHYSICAL REALIZATIONS 

The results presented in this work suggest the intrigu- 
ing possibility of studying the effects of cooperative emis- 
sion with a mesoscopic number of quantum emitters in 
a high-finesse resonator. Moreover, depending on the 
specific physical implementation of the proposed model, 
different regimes of cavity QED might be explored and 
characterized in terms of the statistical properties of the 
emitted radiation. To this purpose, we stress the rele- 



vance of Fig. [3] in providing a roadmap to experimental 
investigations. In fact, depending on the specific system 
of few quantum emitters coupled to a single cavity mode, 
different regimes of n/g might be realized. 

However, some of the requirements that we assumed in 
our analysis might be difficult to be simultaneously real- 
ized in a physical implementation of the proposed model. 
Namely, we assumed identical quantum emitters, equal 
dissipation rates, low pumping of each individual emitter 
as compared to its coupling to the cavity mode. In order 
to realize the conditions 7*, 7 < Px < E, each quan- 
tum emitter should have a high beta factor (quantifying 
the coupling to the cavity loss channel). Usually, this 
is quite challenging with solid state artificial atoms, like 
quantum dots, where inhomogeneous broadening (which 
can be related to 7* ) would inevitably mix the different J 
states. Atomic systems in optical resonators would cer- 
tainly offer an ideal implementation for what concerns 
the atomic-like properties, but it would be difficult to 
have a mesoscopic number of such atoms at fixed spatial 
positions. Alternatively, atomic-like impurities in solid 
state resonators, e.g. semiconductors, might represent a 
valid scenario for this type of physics to be investigated. 

As a first example, we consider semiconductor quan- 
tum dots, especially made in III-V materials, as a promis- 
ing approach to coupling in a deterministic way |22l 155] a 
few quantum emitters to the same cavity mode, e.g. us- 
ing photonic crystal cavities jl4] . Such systems can have 
5 ~ 0.1 meV [221[23| and K ~ 0.01 meV [56]. We can thus 
expect to be in the range K/g ^ 0.1 for such a physical 
implementation, where subradiant vs. lasing behaviors 
could be studied under weak driving conditions. As these 
nanostructures are usually grown randomly, and strong 
inhomogeneous fluctuations between different quantum 
emitters occur, techniques have been recently developed 
to post-process the sample and obtain more homogeneous 
ensembles [57]. Electrical control of different quantum 
emitters has also been shown to provide an effective tun- 
ing mechanism ^58 . Moreover, growth of pyramidal dots 
in pre-patterned positions |59j is a very promising tech- 
nique to place more than one such nanostructures in the 
antinodes of the cavity field in order to maximize cou- 
pling. 

For active centers in solid-state resonators, the cou- 
pling to the cavity mode is usually several orders of mag- 
nitude smaller, dictated by the smaller oscillator strength 
of an atomic-like impurity as compared to a quantum dot. 
In such case, the conditions for high-^/g values might 
be investigated, i.e. the crossover from subradiant to 
superradiant emission, and maybe to individualization. 
Recently, single to few Nitrogen- Vacancy centers in dia- 
mond have been implanted in pre-determined spatial po- 
sitions 6OJ . Resonators in diamond have already been re- 
alized with quality factors ranging from few hundreds in 
early attempts [6T1[62] up to 4000 more recently [63]. Fur- 
thermore, NV-centers are not the only type of active de- 
fects in diamond, see e.g. recent advances with Cr-centers 
[51] and Silicon- vacancies (SiV centers) [jBSl. Among solid 
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state emitters, we also include single molecules in stand- 
ing wave resonators made from a Bragg mirror and a 
resolution limited focused beam through an optical fi- 
bre fBSI . The latter approach has the advantage that the 
cavity can be formed at any arbitrary position, once the 
single molecules have been selected on the Bragg mirror 
surface. Such quantum emitters display radiation-limited 
zero-phonon line, thus showing very limited pure dephas- 
ing rate with almost radiation-limited linewidths p7| . 

Finally, a promising platform to realize a controlled ex- 
periment might be with circuit QED devices, where the 
coupling of a few qubits to a single resonator has already 
been realized [TS] . High tunability of the qubits physical 
properties, such as transition energy, dissipation rates, 
pure dephasing and effective oscillator strength, make 
these systems particularly attractive to closely realize the 
model in Eqs. (l][2), and to experimentally address the 
physical results reported in the present work. 



VII. CONCLUSIONS AND PERSPECTIVES 

We have presented an extensive theoretical analysis of 
the quantum statistical properties of a source of elec- 
tromagnetic radiation made of few two-level emitters 
coupled a single resonator mode. We have assumed a 
model in which the two-level systems are incoherently 
pumped by a weak external drive, and dissipate their en- 
ergy through spontaneous emission in free space, pure 
dephasing, or other loss channels that can be described 
within a Liouvillian formalism. In this sense, the present 
work should be considered as an extension of our previ- 
ous report on single emitter cavity QED [33], where we 
have analyzed the crossover from a single-photon source 
(emitting anti-bunched radiation) to a single-emitter las- 
ing regime. 

However, the case of more than one emitter studied 
here presents a much richer phenomenology depending 
on the specific light-matter coupling regime of each sin- 
gle emitter with the common radiation field, which can 
be identified from the statistical properties of the emit- 
ted radiation. Throughout this work, we have assumed 
a tuning parameter given by the ratio of the cavity de- 
cay rate to the single emitter-cavity coupling strength, 
which we could numerically investigate by spanning over 
10 decades. As a function of this parameter, we have 
clarified the crossover from lasing (low cavity decay rate 
as compared to emitters cavity coupling) to chaotic light 
emission (very large cavity decay rate), going through 
sub- and super-radiance in steady state. Although in the 
superradiant regime the emitted radiation would present 
a close-to-Poissonian statistics, we have explicitly shown 
that it would not reach the level of coherence displayed 
by a laser source. Between these two regimes, sub-radiant 
emission takes place and the emitted radiation undergoes 
super-Poissonian intensity fluctuations. 

We have shown that the strong bunching in the degree 
of the second-order coherence oscillates with the parity of 



the number of emitters, eventually converging to a single 
value in the large number limit. We have clarified the 
origin of this novel behavior by using a suitable basis 
for the system of few emitters, showing that it is due to 
the presence of extra-loops in the radiative cascade from 
subradiant states when the number of emitters is odd, 
thus reducing the degree of bunching with respect the 
even case. Finally, we have also studied the conditions 
for reaching an incoherent source of chaotic radiation, 
which can be done either by coupling the ensemble of 
emitters to a very leaky cavity, or by enhancing the pure 
dephasing rate (a feature that is quite natural for solid- 
state quantum emitters). 

As a last remark, notice that the cavity finesse has 
been the only parameter to be significantly changed in 
order to produce very different output fields. Recently, 
fast dynamical control of this parameter has been demon- 
strated in different scenarios pHl l69] , suggesting that a 
system composed of few emitters inside a resonant cav- 
ity could be a versatile source of very different quantum 
states of light. An immediate and important extension 
of this work related to its experimental and technologi- 
cal implications would be to investigate the effects caused 
by the inhomogeneous distribution of the quantum emit- 
ters, which is a typical setting for smiconductor quantum 
dots. A further seemingly important extension would be 
to explore the consequences of the ultra-strong coupling 
regime (available nowadays for superconducting qubits in 
microwave resonators) on the statistics of the emitted ra- 
diation. In the latter scenario, anti-rotating wave terms 
in the hamiltonian ([T]) must be taken into account [70] in 
order to properly get the emitted radiation and statistics. 
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Appendix A: Basis for the N — 3 case 

We hereby derive the basis states employed for the an- 
alytical estimation of 5^^-'(0) in the three emitters case. 
States with J < Jmax are degenerate as soon as > 2. 
We lift the degeneracy by exploiting the permutation 
symmetry in the free Hamiltonian, Eq. ([T]). Namely, 
we use symmetry-adapted state vectors related to the si- 
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Q] ® [2] ® [3] = 

(j=i) e 

(j=i) e 

(j=f) e 



fip 



A = (2,l) 



A=f2,l) 



1 = 1 



i = 2 



fTT2T3] A = (3,0) 



FIG. 6; Total spin representations of three a 
Young tableaux. 



I spins with 



multaneous diagonalization of the angular momentum J 
(with Jniax — N/2) algebra with the permutation group 
P/v of N atoms [39]. We write the states in the coupled 
basis as 



J A 

Mj i 



(Al) 



where J and Mj represent the collective dipole and its 
projection labels whereas i enumerates the states be- 
longing to the A permutation irreducible representation. 
The derivation reduces to the use of standard branching 
rules for Young tableaux with the theorem in Sec. V of 
Ref. [39] or with the tensor method [40]. Since we are 
dealing with a tensor product of two-dimensional Hilbert 
spaces, ~ C^, then for 3 emitters the branching rules 
give three types of tableaux, out of which we must cancel 
those with more than two rows, as schematically outlined 
in Fig. [6] By means of Young symmetrizers, the coupled 
basis is easily obtained in terms of the uncoupled states 
as 



3 3 
2' 2 
3 1 
2' 2 
3 1 
2'~2 
3 3 
2'~2 

2 



eee) 



{\eeg) + \ege) + \gee)) 



= ~^i\^99) + Igge) + \geg)) 



= 1.935) 



= (2,1) 
= 1 

= (2,1) 
i = 1 



V6 



[2\eeg) - \ege) - \gee)) 



V6 



[^\99e) - \geg) - \egg)) 



2 



(2,1) 

= 2 

= (2,1) 
i = 2 



= y|(l^'5e) - l^ee)) 



V2 



(Iseg) - \egg)) 



Appendix B: Analytic solution for 3 emitters 

In this appendix we show the analytical details for cal- 
culating the autocorrelation at zero-time delay for the 
iV = 3 case. 

The steady state solution to the system of equations 
relevant for the g^^)(0) calculation in the case of 3 emit- 
ters gives: 

P.. 



P3/2,-l/2 
P3/2,l/2 = 
Pl/2,1/2 = 
Pl/2,-1/2 
P3/2,3/2 = 



- -jTP3/2,-3/2 

p,(p, + 2r) 



4r2 

P,(2F, + 5r) 

r(r + 6P,) 
4P^ -t- 3r) 



P3/2,-3/2 
P3/2,-3/2 



(Bl) 



"P3/2,-3/2 



"P3/2,-3/2 ■ 



r + 6P, 

p2(42r2 + 29rP^ + 6P^) 

i2r3(r-f 6P^) 

Once we neglect the off-diagonal terms (see text in 
Sec. IV B[ ), the two degenerate states |l/2, ±1/2; 1), 
|l/2, ±1/2; 2) give the same values, and here they are 
considered simply as 1 1/2, ±1/2). The lowest matrix ele- 
ment is calculated from the constrain on the trace of the 
density matrix. It reads 

12r3(r + 6P2;) 

P3/2,-3/2 



84r4 -I- 306r3p^ + 20ir2p2 + 47rp3 



In powers of P^j/F we finally obtain 



6Pr 
(B2) 



(0) 



14 F 
75 P^ 

-23.4 



2.067- 3. 50F 



P:r. 



9.18 



(B3) 



The same argument could have been used for the two- 
emitters case, for which the solution reads 

2 



4 F 



13 



2 Pr 



0^(0) ~ + — 

^ ^ ' 9P^ 27 27 F 



O 



(B4) 



In the case we are considering, namely T/P^ = 40, 
the error in neglecting all the term s but the first one in 
F/P^ is of about 3% for Eq. ([B4]), while for eq. (|B3]) it 
would have been ~ 30 %. This formally explains why 
in the two-emitters case the linear approximation works 
so well. Indeed, setting V/Px — 40 in Eq. (B3| we have 
^(2)(o) ~ 9.46 as given in the text. 



Appendix C: Individualization and chaotic limit 

Developing the expression for J± with respect to indi- 
vidual atomic operators leads to 



E,,,..,«(0)<(r)a^(r)aL(0)) M 



V 



(CI) 
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We are considering the case where the atomic dipoles 
are independent and incoherently excited. It means that 
whatever i,j we have 

U,) = (6>(0> 
(6) - 0, 

being either <tI_ or o'!|_. It is immediate to see that 
whenever i ^ j, k, I (and permutations) the four op- 
erators expectation value is zero. We are left with 
i — j — k — I and the case where the indices are equal to 
each other in pairs. As a consequence 



a^ = E(4(o)4m^-W'^-(o)) 



contribution as the field emitted by a single atom is an- 
tibunched, the pairs {i = k,j = I) and (i = l,j = k) 
contribute with 2N{N — 1) terms equal to P (by defi- 
nition, we have ^["'^^(O) = 1). Eventually, one finds the 
^(2)(o) = 2(1 - 1/N) result mentioned in the text. It is 
worth noticing that in the large number limit we have 
(7(^^(0) = 2 corresponding to a chaotic chaotic source 

nasi]. 

It is also interesting to remark that with a distribution 
of N classical electromagnetic emitters we still have 



-1 (C2) 



k^l ^ 

+ (a';(0)4(r))(aL(r)aL(0)) 



k k^l 

where Ik = {o'\{Q)cr^{Q)) For zero delay, the denomi- 
nator V involves {J^J_) = ^^^{a\a-'_), which simply 

writes '}2,i{^+'^-) = ^-^1 ^ the dipoles are independent. 
As for the numerator, the case i = j = k ~ I gives no 



However, the classical Schwartz inequality imposes 
Qi (0) > 1 and therefore we now obtain for the N emit- 
ters 

^''^(0)^^+2(1-1)^2-1. (C3) 

The value 2(1 — 1/A^) is therefore somewhere interme- 
diate between the pure particle like behavior, 1 — 1/A^ 
(see main text), and the pure classical undulatory pre- 
diction, 2 — 1/A^ (in agreement with the usual Einstein 
interpretation of the fluctuation formula for the black- 
body radiation). 



[1] R.H. Dicke, Phys. Rev. 93, 99 (1954). 

[2] M. Gross and S. Haroche, Phys. Rep. 93, 301 (1982). 

[3] H.M. Gibbs, Q.H.F. Vrehen, and H.M.J. Hikspoors, 

Phys. Rev. Lett. 39, 547 (1977). 
[4] M. Scheibner, T. Schmidt, L. Worschech, A. Forchel, G. 

Bacher, T. Passow, and D. Hommel, Nat. Physics 3, 106 

(2007). 

[5] R. Bonifacio and L.A. Lugiato, Phys. Rev. A 11, 1507 
(1975). 

[6] J.P. Clemens and H.J. Carmichael, Phys. Rev. A 65, 
023815 (2002). 

[7] R. Bonifacio, P. Schwendimann, F. Haake, Phys. Rev. A 
4, 302 (1971). 

[8] V.V. Temnov and U. Woggon, Phys. Rev. Lett. 95, 
243602 (2005). 

[9] V.V. Temnov and U. Woggon, Opt. Express 17, 5774 

(2009) . 

[10] D. Meiser and M.J. Holland, Phys. Rev. A 81, 033847 

(2010) . 

[11] D. Meiser and M.J. Holland, Phys. Rev. A 81, 063827 
(2010). 

[12] A.N. Poddubny, M.M. Glazov, and N.S. Averkiev, Phys. 

Rev. B 82, 205330 (2010). 
[13] F.P. Laussy, A. Laucht, E. del Valle, J.J. Finley, and J.M. 

Villas-Boas, arXiv:1102.3874 (2011). 
[14] S. Strauf, K. Hennessy, M.T. Rakher, Y.-S. Choi, A. 

Badolato, L.C. Andreani, E.L. Hu, P.M. PetrofT, and D. 



Bouwmeester, Phys. Rev. Lett. 96, 127404 (2006). 
[15] L. DiCarlo, J.M. Chow, J.M. Gambetta, L.S. Bishop, 

B.R. Johnson, D.I. Schuster, J. Majer, A. Blais, L. Frun- 

zio, S.M. Girvin, and R.J. Schoelkopf, Nature (London) 

460, 240 (2009). 
[16] J.P. Reithmaier, G. Sek, A. Loffler, C. Hofmann, S. 

Kuhn, S. Reitzenstein, L.V. Keldysh, V.D. Kulakovskii, 

T.L. Reinecke, and A. Forchel, Nature 432, 197 (2004). 
[17] T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H.M. 

Gibbs, G. Rupper, C. Ell, O.B. Shchekin, and D.G. 

Deppe, Nature 432, 200 (2004). 
[18] A. WaUraff, D. I. Schuster, A. Blais, L. Frunzio, R.S. 

Huang, J. Majer, S. Kumar, S.M. Girvin, and R.J. 

Schoelkopf, Nature 431, 162 (2004). 
[19] E. Peter, P. Senellart, D. Martrou, A. Lemaitre, J. Hours, 

J.-M. Gerard, and J. Bloch, Phys. Rev. Lett. 95, 067401 

(2005) . 

[20] N. Le Thomas, U. Woggon, O. Schops, M.V. Artemyev, 

M. Kazes, and U. Banin, Nano Lett. 6, 557 (2006). 
[21] Y.-S. Park, A.K. Cook, and H. Wang, Nano Lett. 6, 2075 

(2006) . 

[22] K. Hennessy, A. Badolato, M. Winger, D. Gerace, 
M. Atatiire, S. Guide, S. Fait, E. Hu, and A. Imamoglu, 
Nature 445, 896 (2007). 

[23] A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. PetrofT, 
and J. Vuckovic, Nat. Physics 4, 859 (2008). 

[24] M. Kaniber, A. Laucht, A. Neumann, J.M. Villas-Boas, 



12 



M. Bichler, M.C. Amann, and J.J. Finley, Phys. Rev. B 

77, 161303(R) (2008). 
[25] J. Suffczynski, A. Dousse, K. Gauthron, A. Lemaitre, I. 

Sagnes, L. Lanco, J. Bloch, P. Voisin, and P. Senellart, 

Phys. Rev. Lett. 103, 027401 (2009). 
[26] M. Munsch, A. Mosset, A. Auffeves, S. Seidelin, J.-P. 

Poizat, J.-M. Gerard, A. Lemaitre, I. Sagnes, and P. 

Senellart, Phys. Rev. B 80, 115312 (2009). 
[27] J. Kasprzak, S. Reitzenstein, E.A. Muljarov, C. Kistner, 

C. Schneider, M. Strauss, S. Hofiing, A. Forchel, and W. 
Langbein, Nat. Materials 9, 304 (2010). 

[28] F. Haacke, M.I. Kolobov, C. Fabre, E. Giacobino, and S. 

Reynaud, Phys. Rev. Lett. 71, 995 (1993). 
[29] A. Aufleves, J.-M Gerard, and J. P. Poizat, Phys. Rev. A 

79, 053838 (2009). 
[30] C.W. Gardiner and P. Zoller, Quantum Noise (Springer- 

Verlag, Berlin 2004). 
[31] J. Kimble, M. Dagenais, and L. Mandel, Phys. Rev. Lett 

39, 691 (1977). 

[32] J. Dalibard, K. Molmer, and Y. Castin, Phys. Rev. Lett. 
68, 580 (1992). 

[33] H.J. Carmichael, Statistical Methods in Quantum Optics 

vol. 2 ( Springer- Verlag 2008). 
[34] A. Auffeves, D. Gerace, J.-M Gerard, M.F. Santos, L.C. 

Andreani, and J.P. Poizat, Phys. Rev. B 81, 245419 

(2010). 

[35] E. del Valle, F.P. Laussy, and C. Tejedor, Phys. Rev. B 
79, 235326 (2009). 

[36] A. Gonzalez-Tudela, E. del Valle, E. Cancelheri, C. Teje- 
dor, D. Sanvitto, and F.P. Laussy, Opt. Express 18, 7002 
(2010). 

[37] G. Yeoman and G.M. Meyer, Phys. Rev. A 58, 2518 
(1998). 

[38] E. del Valle, F.P. Laussy, F. Troiani, and C. Tejedor, 

Phys. Rev. B 76, 235317 (2007). 
[39] F.T. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, 

Phys. Rev. A 6, 2211 (1972). 
[40] Wu-Ki Tung, Group Theory in Physics (World Scientific 

Publishing Co. Pte. Ltd. 2003). 
[41] Y. Kubo, F.R. Ong, P. Bertet, D. Vion, V. Jacques, D. 

Zheng, A. Drzau, J.-F. Roch, A. Auffeves, F. Jelezko, J. 

Wrachtrup, M. F. Barthe, P. Bergonzo, and D. Esteve, 

Phys. Rev. Lett 105, 140502 (2010). 
[42] L Diniz, S. Portolan, R. Ferreira, J.-M. Gerard, P. Bertet, 

and A. Auffeves, arX iv:1101.1842 [quant-ph]. 
[43] C. Cohen- Tannoudji, Cours du College de France, 1979- 

1980. 

[44] R. Loudon, Rep. Prog. Phys. 43,58 (1980). 

[45] E. Jakeman, E.R. Pike, P.N. Pusey, and J.M. Vaughan, 

J. Phys. A: Math. Gen. 10, L257 (1977). 
[46] H.J. Carmichael, P. Drummond, P. Meystre, and 

D. F. Walls, J. Phys. A: Math. Gen. 11, L121 (1978). 
[47] H. Paul, Rev. Mod. Phys. 54, 1061 (1982). 

[48] R. Brouri, A. Beveratos, J.-P. Poizat, P. Grangier, and 
P. Photon, Opt. Lett. 25, 1294 (2000). 



[49] S.C. Kitson, P. Jonsson, J.G. Rarity, and P.R. Tapser, 

Phys. Rev. A 58, 620 (1998). 
[50] Y.Y. Hui, Y.R. Chang, T.S. Lim, H.Y. Lee, W. Fann, 

and E.G. Chang, Appl. Phys. Lett. 94, 013104 (2009). 
[51] K.D. Weston, M. Dyck, P. Tinnefeld, C. MuUer, 

D.P. Herten, and M. Saner, Anal. Chem. 74, 5342 (2002). 
[52] A. Cuche, O. MoUet, A. Drezet, and S. Huant, Nano Lett. 

10, 4566 (2010). 
[53] S. Pezzagna, B. Naydenov, F. Jelezko, J. Wrachtrup, and 

J. Meijer, New J. Phys. 12, 065017 (2010). 
[54] B.D. Gerardot, S. Strauf, M.J.A. de Dood, A.M. By- 

chkov, A. Badolato, K. Hennessy, E.L. Hu, 

D. Bouwmeester, and P.M. Petroff, Phys. Rev. Lett. 95, 

137403 (2005). 

[55] A. Dousse, L. Lanco, J. Suffczynski, E. Semenova, A. 

Miard, A. Lemaitre, L Sagnes, C. Roblin, J. Bloch, and 

P. Senellart, Phys. Rev. Lett. 101, 267404 (2008). 
[56] S. Combrie, A. De Rossi, Q. V. Tran, and H. Benisty, 

Opt. Lett. 33, 1908-1910 (2008). 
[57] A. Rastelli, A. Ulhaq, S. Kiravittaya, L. Wang, A. Zren- 

ner, and O. G. Schmidt, Appl. Phys. Lett. 90, 73120 

(2007). 

[58] A. Laucht, J.M. Villas-Boas, S. Stobbe, N. Hauke, F. Hof- 

bauer, G. Bohm, P. Lodahl, M.-C. Amann, M. Kaniber, 

and J.J. Finley, Phys. Rev. B 82, 075305 (2010). 
[59] M. Baier, S. Watanabe, E. Pelucchi, E. Kapon, S. Varout- 

sis, M. Gallart, L Robert-Philip, and I. Abram, Appl. 

Phys. Lett. 84, 648 (2004). 
[60] B. Naydenov, R. Kolesov, A. Batalov, J. Meijer, S. Pez- 
zagna, D. Rogalla, F. Jelezko, and J. Wrachtrup, Appl. 

Phys. Lett. 95, 181109 (2009). 
[61] C.F. Wang, Y-S. Choi, J.G. Lee, E.L. Hu, J. Yang, and 

J.E. Butler, Appl. Phys. Lett. 90, 081110 (2007). 
[62] C.F. Wang, R. Hanson, D.D. Awschalom, E.L. Hu, T. 

Feygelson, J. Yang, and J.E. Butler, Appl. Phys. Lett. 

91, 201112 (2007). 
[63] A. Faraon, RE. Barclay, C. Santori, K.M.C. Fu, and R.G. 

Beausoleil, Nat. Photonics, 5, 301 (2011). 
[64] L Aharonovich , S. Castelletto , B.C. Johnson , J.G. Mc- 

Callum, and S. Prawer, New J. Phys. 13, 045015 (2011). 
[65] C. Wang, C. Kurtsiefer, H. Weinfurter, and B. Burchard, 

J. Phys. B 39, 37 (2006). 
[66] C. TonineUi, Y. Delley, T. Stoferle, A. Renn, S. 

Gotzinger, and V. Sandoghdar, Appl. Phys. Lett. 97, 

021107 (2010). 

[67] C. Toninelli, K. Early, J. Bremi, A. Renn, S. Gotzinger, 
and V. Sandoghdar, Opt. Express 18, 6577 (2010). 

[68] P. J. Hardingand, T. G. Euser, Y.-R. Nowicki-Bringuier, 
J.-M. Gerard, and W. L. Vos, App. Phys. Lett. 91, 
111103 (2007). 

[69] J. Upham, Y. Tanaka, T. Asano, and S. Noda, Opt. Ex- 
press 16 21721(2008). 

[70] S. De Liberato, D. Gerace, L Carusotto, and C. Ciuti, 
Phys. Rev. A 80 053810 (2009). 



